Non-Abelian magnetic monopole in a Bose-Einstein condensate 
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Recently, an effective non-Abelian magnetic field with a topology of a monopole was shown to 
emerge from the adiabatic motion of multilevel atoms in spatially varying laser fields [J. Ruseckas 
et al., Phys. Rev. Lett. 95, 010404 (2005)]. We study this monopole in a Bose-Einstein condensate 
(BEC) of degenerate dressed states and find that the topological charge of the pseudospin cancels 
the monopole charge resulting in a vanishing gauge invariant charge. As a function of the laser 
wavelength, different stationary states are classified in terms of their effect to the monopole part of 
the magnetic field and a cross-over to vortex ground state is observed. 
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Introduction. — Existence of an isolated point source 
of magnetic field, that is, a magnetic monopole, has been 
an intriguing question ever since Dirac proposed that the 
existence of a magnetic monopole leads in a natural way 
to the quantization of electric charge 1] . The magnetic 
monopole considered by Dirac is intrinsically a singular 
object with an string-like singularity attached to it. In 
the seminal work of 't Hooft and Polyakov, non-singular 
monopole configurations were shown to arise in the con- 
text of non-Abelian gauge theories [3, y] . Unfortunately, 
these kind of monopoles are expected to be extremely 
heavy which renders their experimental detection dif- 
ficult. Despite ingenious theoretical and experimental 
work in the fields of high energy physics, cosmology, and 
condensed matter physics, convincing evidence of the ex- 
istence of stable magnetic monopoles is still lacking [J| . 

Advances in trapping and manipulating degenerate 
quantum gases have revealed the potential of the cold 
atom systems to serve as quantum simulators for ideas 
beyond the usual condensed matter phenomena. In par- 
ticular, it has been proposed that non-Abelian gauge po- 
tentials are realized in the effective description of atoms 
with degenerate internal degrees of freedom coupled to 
spatially varying laser fields [aja, [Zl, |8| . Based on the ear- 
lier work by Wilczek and Zee Q , Ruseckas et al. Q have 
shown that non-Abelian gauge potentials describe the off- 
diagonal couplings between the degenerate dressed states 
in atoms with multiple degenerate internal states. As a 
specific example, a U{2) gauge potential corresponding to 
a non-Abelian magnetic monopole was constructed ^ . 

Previous monopole studies in Bose-Einstein conden- 
sates (BECs) have concentrated on systems where the 
monopole state occurs without any gauge potential and 
stems essentially from the non-linear interactions be- 
tween atoms in different hyperfine spin states [lO|. This 
renders such monopoles typically energetically unfavor- 
able [lO|. On the other hand, atom-laser interaction 
induced Abelian magnetic monopoles and non-Abelian 
monopoles in optical lattices have been investigated re- 



cently [n|. In this Letter, we consider an isolated non- 
Abelian magnetic monopole in a Bose-Einstein conden- 
sate, which is similar to the monopoles in the Yang-Mills- 
Higgs model [I4]. We study the lowest energy states of 
the system and show that they can be classified accord- 
ing to the topological charge of the pseudospin texture 
generated by gauge transformations. On the other hand, 
we find the gauge invariant total charge of the system to 
be always zero. 

The model. — Consider a four-level tripod scheme, 
in which degenerate atomic states {|1), |2), 13)} are ex- 
cited to a common virtual state |0), see [6|. In the 
rotating wave approximation the atom-light interaction 
Hamiltonian has two degenerate dark states |xi) and 
1X2) corresponding to zero eigenvalue, and two bright 
states Ixg) and \x^ with a non-zero projection to the 
radiatively decaying excited level |0). Expressing the 
full quantum state |<i>) 

the dressed states as |$) — X]m=i '/'m(^)|x 
rise to a unitary transformation ZYmfe(''") = (Xm('")l^) be- 
tween vectors * = (V'l ■ ■ • i'iY' and C = (Co ■ ■ ■ Ca)'^ as 
i^m[r) = X]fc=o^™fc(^)Cfe(^)- If transitions from the de- 
generate dark states to bright states are neglected one 
obtains an effective single-particle Hamiltonian for the 
dark states 



^k=o^k{r)\k) in the basis of 
^'- '-^'^ (r)) gives 



Hi= dr 



— {D^^j,)\D,xj,) + t/>t(Fe.t + $)^ 



, (1) 



where D^^ = d^ — iA^^ and ip — (V'i,V'2)"^, see Ref. [y]. 
The Berry connection A ~ A^e^ and the scalar potential 
$ transform as A -> UAU^ - i{VU)W and $ -^ U^W 
under a local change of basis given by a unitary matrix 
U{r). Furthermore, we denote A^ — A^cr" and a°' — 
{l,a^,a''^ ,a^) [19|. The external potential that confines 
the atoms is given by Voxt- 

The non-Abelian magnetic monopole was shown to 
arise from two circularly polarized beams and a linearly 
polarized beam yielding [6| 



mean field Hamiltonian 



cosz9 ^ 1,, ,, , 



n= dr 



X [(1 + 008^^)1 + sin^iJcr^], (2) 



$ = 
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-irl- -—sin ip sin 2i? ct^ 
r 2r 



fc2 + /c'- 



sin2 2t?(l-cr^) 



(3) 



where (r, (p, t?) denote the spherical coordinates. As 
pointed out in [g, the effective magnetic field B^ = 
\^ ^Lv\^u\^^ is of the form B — 1/r^ e^-a^ + • • • in 
which the omitted "non-nionopole" terms do not gen- 
erate a net flux through a closed surface enclosing the 
monopole. The field strength tensor F" is given by 
F^u = d^Al^d,Al+2e'^^'A'^^Al andF°, ^ d^Al-d^A^, 
see [19]. 

Let us consider the atom-laser and the atom-atom in- 
teractions in the original basis {|A;)}, fc = 0, • • • ,3. For 
a system consisting of, e.g., *^Rb atoms with hyperfine 
spin F, the set of degenerate internal states {|1), |2), |3)} 
can be taken to be the Zeeman sublevels \F — l,mp), 
rrip — 1,-1,0 at the ^St^i^ multiplet, and the excited 
state e.g. \F = 0,mp — 0) at the 5P3/2 multiplet. For 
®^Rb the scattering lengths in the channels corresponding 
to total hyperfine spin f — and f = 2 of the two scatter- 
ing particles are almost equal [l3|, and hence we neglect 
the spin-dependent part of interaction between states 
\F = l,77iy), see [IJ]. Furthermore, similarly to Ref. |8[ 
we assume that the population of the two bright states 
Ixa) and Ixi) is vanishingly small compared to the popu- 
lation of the dark states \xi) and \x2)- This is, in fact, a 
non-trivial assumption since for a three-level A-system it 
has been shown that the dark state can be unstable un- 
der small deviations, leading to non- vanishing population 
of the bright state [l5(. Similar analysis for the tripod 
system considered here has yet to be performed. Using 
the unitary transformation ZYmfe(r) and the Rabi frequen- 
cies i ni,2{r) = rio^e'C^^T^), ft^^r) = r!o|e'*'-'" cor- 
responding to the monopole field, one obtains an effective 
Hamiltonian 7?^^ — ^ J dr {ij}'^ ij})'^ describing the inter- 
actions in the dark state manifold. 

In the absence of light induced gauge potentials Bose- 
Einstein condensation takes place at low enough tem- 
peratures. On the other hand, the existence of artificial 
gauge potentials can alter the formation of BEG, leading 
to fragmentation [l6|. Although strictly speaking the 
mean field approach used below cannot treat effects such 
as fragmentation, even in the absence of a true conden- 
sate the mean field states can be used to construct ap- 
proximations for the true ground state as well as to inter- 
pret the possible outcomes of a particular experimental 
realization 17[. Replacing xjj by {tj}) = ip we obtain a 



2m 



{D^xP)\D^'iP)+xP\V,^, + <P)xP+'^\'iP\' 



(4) 

The atom-laser interaction is diagonal in the basis con- 
sisting of the dressed states {|x('"))}- The dark states 
correspond to a zero eigenvalue and hence the atom-laser 
interaction does not appear in Eq. fl]). For optically con- 
fined spin-1 bosons such as ^^Rb, the external trapping 
potential Vext can be treated as a scalar function [1J| and 
the model in Eq. ([4]) is U{2) gauge invariant. 

Gauge invariance. — In the Abelian case, the magnetic 
monopole is accompanied with singular filament, the 
Dirac string, that extends outwards from the monopole. 
For non- Abelian monopoles it has been pointed out that 
the Dirac string can be removed by a suitable gauge 
transformation 12]. Similar line of reasoning applies also 
here and by applying a gauge transformation Uo{r) — 
^iipcr ^-i-da'" II Qj^g gg^j^ remove the divergent part of A. 
The transformed vector potential is given by 

A! = -eh'^a" - ^eJcos(2ip)a« ~ sin(2(^)CT^] 
r ^ 2r 

+ Uke, - fc'e^)[(l + cos^ d)l + sin^ d feV], (5) 

with h = (sin I?, sin(2(^) cosi?, cos(2iy9) cosi?). The scalar 
potential transforms according to $' = {/o$[/q . It should 
be noted that A' is not well-defined at the z axis al- 
though the gauge transformation removed the explicit 
divergence. The scalar potential <&' on the other hand, 
remains well-defined at the z axis. 

The topological structure of the condensate can be an- 
alyzed by studying the pseudospin 



s{r) = xp* {r)cnp{r). 



(6) 



We note that s is not itself gauge invariant, but trans- 
forms under any gauge transformation [/ as s° — > u'^^s^ 
where U'^a°'U = u°'^a^ . An explicit computation of the 
matrix elements m°^ shows that the gauge transformation 
of the pseudospin s" can be written as 

s'^ct" -^ Us^a^U^. (7) 

We define the covariant derivative of the pseudospin by 

P^s'^ = iP^a-'D^iP + c.c. = d^s"- + 2£'^''=4s^ (8) 

from which we see that P^s" transforms as P^s" — > 



M^^P^S^ 



Using identities 2|s| 



trfs'^cr'^sV'') 



2gafccgap^g62?^^c ^ -i tr[s°(T° (X'^s'')(T''(X>^s'=)cr'=] , and 
Eq. (O one observes that 

1 

isl^-" 2^ 






e''''=s°I?^s^I?,s^ (9) 



is invariant under U{2) gauge transformations of xp and 
Afj_. Tensor G^j/ is the "electromagnetic" tensor intro- 
duced by 't Hooft in the studies of monopoles in unified 



gauge theories 21 , although here s is not the fundamental 
variable but defined by Eq. ([6]). 

Using a unit vector field n = •s/|s| one can write G^^ 



in the form [l2] G^^ 



M„ 



Hu,p where 



M^^^d^in-^AD^d^n'^Al), (10) 

1 ....„„ ... .. ^^^^ 



^A- - ^e-''"'h''d^,n''d,n 



In particular, Eqs. pII)) - pT|) imply that G^^, is linear 
with respect to the gauge potential AJ^. This enables 
us to study separately the terms in the gauge potential 
which are responsible for the monopole and consider the 
terms depending on k or k' as a background field for the 
monopole configuration, see Eq. ([5]). Using the tensor 
G^u we define magnetic charge density J — ^n,j\d^i,Gu\, 
and with the Gauss' theorem the total magnetic charge 
can be written as [I4I 

^^^ I '^"■^ ^M-A [Mux - H^x] - Qm - Qs- (12) 

where dr^ is the surface element on 5°°. The factor 
of 2 in the derivative of s" [Eq. ([8])] implies that Qg is 
quantized in the units of i, that is, Qg is i times the 
winding number of the unit vector field n, see [12l |. 
The gauge transformation Uq can be decomposed as 



Uo 
Ui 



U2U1 where and U2 



^iva^ ^Mayi2^-%Va^ ^^^ 



^iipa 



In particular, transformation with Ui brings 
the monopole part of A to the form which has the Dirac 
string at negative z axis and is thus analogous to the 
Dirac description of the Abelian monopole il!|. The 
meaning of the different gauges can be understood qual- 
itatively by considering the following ansatz '0±(r) = 
<j>{r)C± with C± = (1/a/2, T1/a/2)^ for the condensate 
order parameter in the non-singular gauge Q2 = {A' , $'). 
For this simple trial wavefunction one obtains the follow- 
ing diagram 



Q'ii^^o 



Qf 



)(i) 



= ±1 







Qf^±l 



Go 



Q 



(0) _ 



(13) 



0, 



where Qq = {A, $) is the original gauge given by Eqs. ^ 
and ([3]) and Qi is obtained from Qo with Ui. In the next 
section we will use the intermediate gauge Gi to classify 
different stationary states. 

Numerical analysis. — To study the lowest energy 
states we search numerically solutions to the field equa- 
tions corresponding to the Hamiltonian in Eq. ^. We 
discretize the Hamiltonian using the scheme previously 
utilized in the lattice gau ge theory studies of the non- 
Abelian Higgs model [l8|. We replace {D ^ipy (D ^xp) 

by G^,x = {rpx+f. - U^^x'^xVi'fPx+f, " f^^.x'^x)/^', 
where X — {u,v,w) G N"^ denotes one lattice site, 
U, 



constant a tends to zero, Cp,,x reduces to {D f^ipy [D ^xp) . 
The rest of the Hamiltonian is readily discretized yielding 
the mean field energy 



^-^ = £ E ^P,^ + E[^i(^-t,x + 1>x)V'x + y IV'xh 



X,p 



X 



Minimizing the free energy Td ^ £d - Mch Z^jc I'^xP 
gives the discrete version of the time- independent Gross- 
Pitaevskii (GP) equation 



E e^Ix^x+A - S-^x + U^,x-[.'*l'x-ii. 



2ma^ 



+ (14xt,x + ^x)'>Px + col-f/'xl Vx = Mch-^x- (14) 

The chemical potential guaranteeing the conservation of 
the particle number is denoted by fi^^i- 

For the numerical calculation we assume the external 
potential of the form l^xt = ma;^r^/2 and for simplicity 
that the wave vectors of the laser fields are equal and 
A: = fc' > 0. The spatial variables are scaled with the 
harmonic oscillator length a^^^ — ^Jh/muj^. The strength 
of the inter-particle interaction is taken to be c^N /huj^ — 
100. We solve the GP equation P^ iteratively using the 
successive over-relaxation scheme (SOR) with periodic 
boundary conditions in a grid of 121^^ points. Since the 
vector potential in Eq. ([5]) is not well defined at the z axis 
we shift the grid in x-direction by a positive constant 
77 taken to be 77 = a/10. The results are independent 
of the choice of 77. The mean field energy of different 
stationary states is shown in Fig. [T] as a function of ka^^^. 
The calculations are carried out in the non-singular gauge 
Q2 ~ (A', $') for which the corresponding charges are 
denoted by Qg and Q\i ■ 



l(2) 



All stationary states have Qg — Q\.j = for fca,^^ < 
3.5. For A:a,j^ > 3.5 the ground state contains a vortex 
and the classification in terms of Eq. (|12p is no longer 
valid. For ka^^^ < 0.055 there is a pair of solutions which 
arc degenerate within the numerical accuracy. These so- 
lutions behave under the gauge transformation C/J such 
-^ Qg = ±1, where the charges in the 



that Q^g'' 



(1) 



i(i) 



^iaA',^ 



intermediate gauge Gi are denoted by Qg ' and Q\l ■ For 
^^ho ^2.1 the pair appears again as excited states with 
Qg — 1 state having slightly larger energy. The differ- 
ence in the energies between these two states is, however, 
very small and cannot be easily distinguished from Fig. 1. 
In the original gauge Go, we find Qg — for all found 
states, similar to Eq. P^ . For ka-^^ < 3.5 the lowest 
energy state is always the state with Qg = 0. 

In the gauge Go, the single-particle operator /iq of the 
mean field Hamiltonian H — J dr [ip^'ho'^ + ^|'0|'*] has 

the symmetry h(){ip) — a^hQ{—(p)a^ . The two-fold sym- 

(1) 



' X — e*" ^".-^ is the link variable, and p. is one lattice metry explains the degeneracy of Qg = ±1 states for 
unit in the direction fi in the lattice. When the lattice small ka^^^ and suggests that the small energy difference 




(c)^^ (d) 



FIG. 1; (Color online) Left panel: 
mean field energies corresponding to 
the solutions of the Gross-Pitaevskii 
equation. States with Qg = ±1 are 
practically indistinguishable in the 
plot. Right panel: isosurfaces of the 
particle densities for fca^o = 0.86 (a), 
fca,„ = 2.6 (b), fea,„ = 5.0 (c), and 
fca,„ = 6.3 (d). The correspond- 
ing points are marked with • in the 
energy diagram. In panel (d) the 
ground state has a vortex. The parti- 
cle density increases from dark (blue) 
to light (green) color. 



between these two states for large ka^^^ can be due to the 
discrete lattice that breaks the symnietry with respect to 
rotations about the z axis. For small fcoj^^, the energies of 
Qg = ±1 and Qg — states almost coincide and they 
are indistinguishable in Fig. 1. The transition from the 
ground state with Qg = to the vortex ground state 
occurs when initially far apart separated pair of coreless 
vortices starts to become more and more tightly bound 
with increasing ka^^^ eventually forming a singular vor- 
tex. Once the vortex state is formed, its energy seems to 
depend only weakly on ka^^^. 

Within the mean field theory employed here, we cannot 
say anything definite about the possible fragmentation 
of the condensate. If the possible fragmentation is in 
this case related to restoring the broken symmetry 17|, 
the true ground state could correspond to a fragmented 
condensate due to the two-fold symmetry of the mean 
field Hamiltonian. Numerically, however, we seem to find 
only a non-degenerate ground state. To study whether 
the BEG in a monopole field is fragmented, one needs to 
perform a calculation analogous to the one in Ref. [iq . 

Assuming that the underlying system consists of ^^Rb 
atoms, the two obvious gauge invariant observables are 
the particle density g and the hyperfine spin of the con- 
stituent atoms. Both of these quantities can be imaged 
accurately with the state-of-the-art techniques. In the 
original gauge Go, the pseudospin s can be expressed in 
term these two quantities which in principle enables re- 
construction of s. Thus it should be experimentally pos- 



i(i) 



±1 and 



sible to distinguish between states with Qg' 
Qg = by performing a local change of basis using 
Ui. In the current experiments with *^Rb, one typically 
has 1.5 < fcfljj^ < 10 for which the ground state either 
has Qg -^ = or contains a vortex. For very tight traps 
with ujj, of order of kilohertz, the point where Q^ = ±1 

and Qg = states start to coexist can be reached. By 
changing the trap parameters it is then possible to tune 
the value of ka^^^ and explore different regions of the mean 
field phase diagram in Fig. 1. In particular, changing lu^ 
non-adiabatically could result in the Qg — ±1 state. 



In conclusion, we have studied a non-Abelian monopole 
in a Bosc-Einstein condensate and showed that the sys- 
tem can be described with an effective U{2) gauge invari- 
ant model. We identified a gauge invariant charge char- 
acterizing the system and classified different stationary 
states using this charge. Numerical calculations showed 
that the existence of a monopole in the non-Abelian 
gauge potential gives rise to a pseudospin texture with a 
topological charge that cancels the monopole charge. 

The authors would like to thank M. Nakahara for dis- 
cussions and Jenny and Antti Wihuri Foundation, Emil 
Aaltonen Foundation, and the Academy of Finland for fi- 
nancial support. CSC, the Finnish IT center for science, 
is acknowledged for computational resources. 



[11 


P. 


[21 


G 


[31 


A 


[41 


K 


[51 


K 


[61 


J. 


[71 


L. 


[81 


X 


[91 


F. 


[101 


H 




J.- 



[Ill 

[12] 
[13] 

[14] 
[15] 

[16] 
[17] 



[18] 



A. M. Dirac, Proc. Roy. Soc. Lond. A 113, 60 (1931). 

't Hooft, Nucl. Phys. B 79, 276 (1974). 

M. Polyakov, JETP Lett 20, 194 (1974). 

A. Mihon, Rep. Prog. Phys. 69, 1637 (2006). 

Osterloh et al, Phys. Rev. Lett. 95, 010403 (2005). 

Ruseckas et al, Phys. Rev. Lett. 95, 010404 (2005). 
■H. Lu and Y.-Q. Li, Phys. Rev. A 76, 023410 (2007). 
-J. Liu et al, Phys. Rev. Lett 98, 026602 (2007). 

Wilczek and A. Zee, Phys. Rev. Lett 52, 2111 (1984). 

T. C. Stoof et al, Phys. Rev. Lett 87, 120407 (2001); 

P. Martikainen et al, Phys. Rev. Lett 88, 090404 
(2002); C. M. Savage and J. Ruostekoski, Phys. Rev. A 
68, 043604 (2003). 

P. Zhang et al, Eur. Phys. J. D 36, 229 (2005); J. R. 
Hoh and J. M. F. Gunn (2007), eprint: arXiv:0707.0985. 
J. Arafune et al, J. Math. Phys 16, 433 (1975). 
E. G. M. van Kempen et al, Phys. Rev. Lett. 88, 093201 
(2002). 

T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 
A. P. Itin and S. Watanabe, Phys. Rev. Lett. 99, 223903 
(2007); A. P. Itin et al, Phys. Rev. A 77, 043610 (2008) 
T. D. Stanescu et al, Phys. Rev. A 78, 023616 (2008). 
Y. Castin and C. Herzog, C. R. Acad. Sci. Paris, Serie 
IV 2, 419 (2001); E. J. Mueller et al, Phys. Rev. A 74, 
033612 (2006). 

E. Fradkin and S. H. Shenker, Phys. Rev. D 19, 3682 
(1979). 



[19] The Greek indices ^,v... refer to tlie spatial coordinates peated indices is implied, 

and Q runs over values 0, • ■ ■ 3. Indices a,b... refer to 
tiie Pauli matrices (t — {a^ ,a^ ,a^). Summation over re- 



